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Abstract 

(f) ■ In the framework of perturbative QCD approach, we calculate the branching ratio and CP 

> ; 

(NJ \ asymmetry for B®(B S ) — > -k^K^ and B S (B S ) — > n K°(K ) decays. Besides the usual fac- 



X 



torizable diagrams, both non-factorizable and annihilation type contributions are taken into 



account. We find that (a) the branching ratio of B®(B S ) — > ir^K^ is about (6 — 10) x 10 6 ; 

in ■ 

Br(B s (B s ) — ► n K (K )) about (1 — 3) x 10 7 ; and (b) there are large CP asymmetries in the 



two processes, which can be tested in the near future LHC-b experiments at CERN and BTeV 



45 

experiments at Fermilab. 
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I. INTRODUCTION 



The rare charmless B meson decays arouse more and more interest, since it is a good 
place for testing the Standard Model (SM), studying CP violation and looking for possible 
new physics beyond the SM. Since 1999, the B factories in KEK and SLAC collect more 
and more data sample of rare B decays. In the future CERN Large Hadron Collider 
beauty experiments (LHC-b), the heavier B s and B c mesons can also be produced. With 
the bright hope in LHC-b experiments and BTeV experiments at Fermilab, following a 
previous study of B s — > 7r + 7r~ decay [lj], we continue to investigate other B s rare decays. 

The most difficult problem in theoretical calculation of non-leptonic B decays is the 
calculation of hadronic matrix element. The widely used method is the factorization 
_h (FA) Q. It is a great success in e X p.aini„ g the branching ratio of m any decays ft 
|4|, although it is a very simple method. In order to improve the theoretical precision, QCD 
factorization |5| and perturbative QCD approach (PQCD) Q are developed. Perturbative 
QCD factorization theorem for exclusive heavy-meson decays has been proved some time 
ago, and applied to semi-leptonic B — > D{ir)lv decays 6], the non-leptonic B — > Rn 0, 



7T7T |8[ decays. PQCD is a method to factorize hard components from a QCD process, 
which can be treated by perturbation theory. Non-perturbative parts are organized in 
the form of universal hadron light cone wave functions, which can be extracted from 
experiments or constrained by lattice calculations and QCD sum rules. More information 
about PQCD approach can be found in 

In this paper, we would like to study the B Q a {B a ) -> i^K* and B S (B S ) -> ir°K°(K ) 
decays in the perturbative QCD approach. In our calculation, we ignore the soft final 
state interaction because there are not many resonances near the energy region of B s 
mass. Our theoretical formulas for the decay B s — > nK in PQCD framework are given 
in the next section. In section ILTTj we give the numerical results of the branching ratio of 
B s — > 7iK and discussions for CP asymmetries and the form factor of B s — > K etc. At 
last, we give a short summary in section ITVl 



II. PERTURBATIVE CALCULATIONS 





H cS = ^ <j VudVub [CMOriii) + CMOM] - VtlVt^Ci^O^) j> , (1) 



For decay B s — > nK, the related effective Hamiltonian is given by 
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where Cj(ju)(i = 1, ■ • • ,10) are Wilson coefficients at the renormalization scale fi and 
Oi(i = 1, ■ ■ ■ , 10) are the four quark operators 

Oi = (biU^v-Aiujd^v-A, 2 = {biUi) V - A {.u j d j ) v _ A , 

3 = (bidi) V - A Y^ q {Wj)v-A, A = (Mj)v-A Y, q (Wi)v-A, 

5 = {bidi)v-AY^q{Wj)v+A, 6 = (bidj) V -AY,q(qj<li)v+A, ( 2 ) 

7 = l{ptd l ) V -AY. q e q(.Wi)v+A, 8 = \{pid j )v-AY, q e qi.mi)v+A, 
9 = \(bid i )v- A Y. q e q(.Wi)v-A, 0io = Ubidj)v-AJ2 q e q(qjqi)v-A- 
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Here i and j are SU(3) color indices; the sum over q runs over the quark fields that are 
active at the scale \x = 0(rrib), i.e., q G {u, d, s,c, b}. Operators Oi,0 2 come from tree 
level interaction, while 3 , 4 , 5 , 6 are QCD-Penguins operators and 7 , 8 , 9 , O w 
come from electroweak-penguins. 

Working at the rest frame of B s meson, we take kaon and pion masses Mk ~ M n ~ 0, 
which are much smaller than M Bs . In the light-cone coordinates, the momenta of the B s , 
K and 7r can be written as : 

Mr Mr Mr 

Pi = ^|(l,l,0 T ), P 2 = -|(0,1,0 T ), P 3 = ^|(l,0,O r ). (3) 

Denoting the light (anti-)quark momenta in B, K and 7r as k±, A; 2 and ks, respectively, 
we can choose: 

h = {xxpl, 0, k 1T ), k 2 = (0, X2P2, k 2T ), h = {x 3 pt, 0, k 3T ). (4) 

In the following, we start to compute the decay amplitudes of B s — > nK. 

According to effective Hamiltonian we draw the lowest order diagrams of B s — > 71 K 
in Fig. ^ Let us first look at the usual factorizable diagrams (a) and (b). they can give the 
B s — > K form factor if take away the Wilson coefficients. The operators 0\, O2, O3, O4, Og 
and O10 are (V — A)(V — A) currents, and the sum of their contributions is given by 

I* 1 poo 

F e [C] = 167rCVM| / dx\dx2 / 61^6162^2 
Jo Jo 

x { [(2 - x 2 )0^(x 2 ) -r K (l- 2x 2 )<p P K (x2) 
+ r K (l - 2x 2 )(/) T K (x 2 )]a s (t 1 a )h a (xi, 1 - x 2 , h, b 2 ) exp[-S B (tl) - S K (tl)]C(tl) 

+ 2r K (j) P K {x 2 )a s {tl)h a {l - x 2 , x h b 2 , h) exp[-S B (t 2 a ) - S K (tl)]C(t 2 a )}, (5) 

where r n = m 07T /m B = ml/[m B (m u + m d )], r K = m 0K /m B = m 2 K /[m B (m s + m u )\. 
Cf = 4/3 is the group factor of the SU(3) C gauge group. The expressions of the meson 
distribution amplitudes 4>m, the Sudakov factor Sx(ti)(X = B s ,K,tt), and the functions 
h a are given in the appendix. In above formula, the Wilson coefficients C(t) of the 
corresponding operators are process dependent. 
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The operator 5 , 6 , 7 , 0$ have the structure of (V — A)(V + A), their amplitude is 

n 1 /*oo 

F[[C] = 32irC F M B r 7T / dx x dx 2 / b 1 db 1 b 2 db 2 <Ma?iA) 

x {[0^(^2)-^(x 2 -3)0^(a: 2 ) 
+ Ml - x 2 )^(x 2 )]a a (0^«(a:i, 1 " *2, h, b 2 ) exp[-S B (tl) - S {t\) 

+ 2r K ^(x 2 )a s (t 2 a )h a (l - x 2 , x u b 2 , h) ex P [-S B (t 2 a ) - S K {t 2 a )]C{t 2 a )} . (6) 

For the non-factorizable diagrams (c) and (d), all three meson wave functions are 
involved. Using 8 function 5(bi — b 3 ), the integration of bi can be preformed easily. For 
the (V — A)(V — A) operators the result is: 

32 f 1 f°° 

M e[C] = -—-kC f ^/2N c M 2 b / d Xl dx 2 dx 3 / b 2 db 2 b 3 db 3 4> B {xiM) 
3 Jo Jo 

x { [(x 3 - l)tf(x 3 )<t>£(x 2 ) + r K (l - x 2 )^{x 3 )ct> p K {x 2 ) +r K (l- x 2 )4t(x 3 )<&(x 2 )]C(tl) 
a s (t 1 c )h^(x 1 ,x 2 ,x 3 ,b 2 ,b 3 )exp[-S B (t 1 c ) - S n {tl) - S K {tl)} - [{x 2 -x 3 - l)tf(x 3 )(f)£(x 2 ) 
+ r K (l - x 2 )^(x 3 )^(x 2 ) - r K (l - x 2 )^(x 3 )(j) T K (x 2 )}C(t 2 c ) 

a s {t 2 c )h¥\x u x 2 ,x 3 MM)zM-S B (t 2 c ) - S^tl) - S K (t 2 c )}}. (7) 

For the(V — A)(V + A) operators, the formula is: 

M e[C] = -—nC F >/2N' c M B r n / dx x dx 2 dx 3 \ b 2 db 2 b 3 db 3 (j> B {xi,b 3 ) 
3 Jo Jo 

x { [r K (x 2 +x 3 - 2)(f)^(x 3 )(j)^(x 2 ) - r K (x 2 - x^^iyX^^x^ - r K (x 2 - x 3 )(f)l(x 3 )(j)^(x 2 ) 
- r K (2 -x 2 - x 3 )<t>T(x 3 )<f> T K (x 2 ) - (1 - x 3 )<f>r(x 3 )<f>£(x 2 ) - (1 - x 3 )<ff(x 3 )<i>i(x 2 )]C(tl) 
a s (t 1 c )h^ ) (x 1 ,x 2 , x 3 , 62, 63) exp[-S B (tl) - S v (t].) - S K (t].)] + [r K (l -x 2 + x 3 )0f (xg^j^xa) 
+^x(a;2+^3-l)0^(^3)0x(^2)-rx(a;2+a;3-l)0l , (^3)0x(^2)-^(l-x 2 +X3)0^(x3)0^(x 2 ) 

+ x 3 ^(x 3 )(j)i(x 2 ) - x 3( fi(x 3 )<t)£{x 2 )}C(t 2 c ) 
a s (tl)h?\ Xl ,x 2 ,x 3 ,b 2 ,b 3 )e W {-S B (t 2 c ) - S„(t 2 c ) - S K (t 2 c )}}. (8) 

Similar to (c),(d), the annihilation diagrams (e) and (f) also involve all three meson wave 
functions. Here we have two kinds of amplitudes, M a is the contribution containing the 
operator of type (V — A)(V — A), and is the contribution containing the operator of 
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type (V - A)(V + A). 

M a[C] = -—nC F ^/2N~ c M 2 B / dx x dx 2 dx 3 \ M&i b 2 db 2 4>b{xiM) 

Jo Jo 

x { [ x z<t>i (x 3 ) 4>k{x2) + V K (2 + x 2 + x 3 )(j)^(x 3 )(j)^(x 2 ) - r T r K (x 2 - x 3 )(f) P {x 3 )(t) T K {x 2 ) 
- r\r K (x 2 - x 3 )(j)l(x 3 )(j)K(x 2 ) - r\r K {2 - x 2 - x 3 )(j%(x 3 )(j)K(x 2 )]C(tl) 
a s (t 1 e )h < ^\x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) exp[-5' B (^) - S„(tl) - S K (tl)] - [x2^{xz)^k{x 2 ) 
+ r n r K (x 2 + x 3 ) <p P (x 3 ) <px (x 2 ) + r n r K (x 2 - x 3 ) <p p (x 3 ) 0^ (rr 2 ) + r n r K (x 2 - x 3 ) (pi (x 3 ) 0^ (x 2 ) 

+ r 7T r K (x 2 + x 3 )<i)l(x 3 )(i)K(x 2 )]C(t 2 e ) 
a s {t 2 e )hf\x u x 2 , x 3 , b u b 2 ) exp[-S B (t 2 e ) - S,(t 2 e ) ~ S K (t 2 e )]}, (9) 

M a[C] = -—nC F ^/2N~ c M B \ dx x dx 2 dx 3 \ b x db x b 2 db 2 4> B {xiM) 

Jo Jo 

x {[r K (2 - x 2 )(j)^(x 3 )(f) P (x 2 ) + r K (2 - x 2 )(j)^(x 3 )(j)l(x 2 ) - r n (2 - x 3 )(j) P (x 3 )(j)i(x 2 ) 

-r w (2-x 3 )<f>Z(x 3 )<f>A(x 2 )]C(tl) 
a s (t 1 e )h^(x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S B (t 1 e ) - 5 7r (^) - S K (tl)\ + [r K x 2 ^{x 3 )(t) P K {x 2 ) 
+ r K x 2 <p^(x 3 )(t) T K (x 2 ) - r\x 3 <p p (x 3 )(t)j < (x 2 ) - r n x 3 <pl(x 3 )<pK(x 2 )}C{t 2 e ) 

a s (t 2 e )hf\x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S B (t 2 e ) - S„(t 2 e ) - S K (t 2 e )}}. (10) 

The factorizable annihilation diagrams (g) and (h) involve only two light mesons wave 
functions. F a is for (V — A)(V — A) type operators, and F p is for (V — A)(V + A) type 
operators: 

F a [C] = IQitC f M b / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x { [-x 2 (f)^(x 3 )(f)^(x 2 ) - 2r n r K (l + x 2 )(j) P (x 3 )(l)^(x 2 ) + 2r n r K {\ - x 2 )^{x 3 )(t) T K {x 2 )\ 
^s{t 1 g )h g {x 2) x 3 , b 2 , b 3 ) exp[-S v (tl) - SA-(tJ)]C(*J) 
+ [x 3 <f)^{x 3 )(f)i{x 2 ) + 2r n r K {l + x 3 )4> p (x 3 )(f)^(x 2 ) - 2r n r K {l - x 3 )4>l{x 3 )(t> P K {x 2 )] 

C(t 2 g )a s (t 2 g )h g (x 3 , x 2 , b 3 , h) expi-S^tl) - S K (t 2 g )]}, (11) 
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pi poo 

F[[C] = 32ttC f M| / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x {[r K x 2 (f)^{x 3 )(t)^{x2) - r K x 2 <p^(x 3 )<p T K (x 2 ) + 2r n <p p (x 3 )<pj c (x 2 )} 

a s (t 1 g )h g (x 2 , x 3 , b 2 , b 3 ) exp[-S v (fy - S K (typify 
+ [2r K (p^(x 3 )(f)K(x 2 ) +r 7T x 3 (p^(x 3 )(pK(x 2 ) - r^x 3 0^(x 3 )</4(x 2 )] 

C(^)« s (^)^(^3,a;2,&3,&2)exp[-^(^) - (12) 

From Equation ff3|)- ()12|) . the total decay amplitude for B s — > n + K~ can be written as 



A{Bl 



7T 



~K~ 



= UF t 

-Uv; b v td F e p 



V ud V: b {\c x + C 2 ) - V*V td (~C 3 + C 4 + l -C 9 + Cio) 



~Cg + C*6 + -C7 + C*8 



+ M e [KdV^Ci - ^F td (C 3 + Cg)] 



(( c 7 ) i ;„.u, ( c 8 - ic 9 ) - vsv«Mf (c 5 - ^c 7 



-fBV t * b VtdF a 



1 11 

— C 3 C4 — — Co — — Cm 

3 6 2 



f B V t lV td F a p 



— Cc; -|- Cr — ~Cl — _ Ca 

3 6 2 



and the decay width is expressed as 

T(B° S -> 7T+A- 



128tt 1 v 3 



7T 



(13) 



(14) 



The Wilson coefficient C's should be calculated at the appropriate scale t which can be 

fl 

found in the Appendix of Ref. |8[. The decay amplitude of the charge conjugate channel 
B° s — > n~K + can be obtained by replacing V ud V* b to V* d V ub and V t lV td to V tb V t * d in Eq. lfTBl . 
For the decay B s — > 7r°A , its amplitude can be written as 



A(B° S 



7T 



°K°) 



v ud v: b {c x + l -c 2 ) - v: b v u (- l -c 3 - c 4 + l -c 9 + ic M ) 



— Ch — Ck -|- — C7 -|- — c« 
3 6 2 8 



vwk^a - v*v td (-c 3 + -c 9 ) 



-y; fe ^Mf ( -c 7 - c 5 



V*V td M a [ -c 9 - c 3 



v: b v td Mi ( -c 7 - c f 



—fBV t * b V td F a 



1 11 

— — C4 + -Cg + -Cio 

O O 2 



f B v: b v td Fi 



— 77C5 — C*6 + 77C7 + -Cg 
o 2 



(15) 
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and the decay width is then expressed as 



r 2 /\/f 3 

T(B° a - vr ^) = ^ \A(B° -> vr^°)| 2 . (16) 



III. NUMERICAL EVALUATION 

The following parameters have been used in our numerical calculation jlll . Il2j | : 

M Bs = 5.37 GeV,m ^ = 1.4 GeV,m oi ^ = 1.6GeV, Agf c = 0.25 GeV, / Bs = 230 MeV, 

U = 130 MeV, f K = 160 MeV, t b o = 1.46 x 10- 12 s, \V t * b V td \ = 0.0074, = 0.0031. 

(17) 

We leave the CKM phase angle a = <p2 as a free parameter, whose definition is 

« = arg -77^777 • ( 18 ) 

In this language, the decay amplitude of B s — > ir + K~ in eq.(|13|> can be parameterized as 

A = V: b V ud T - V t * b V td P = V* b V ud T[l + ze« a+ % (19) 

where z = \V t * b V td /V* b V ud \\P/T\, and 5 is the relative strong phase between tree diagrams 
T and penguin diagrams P. z and 5 can be calculated from PQCD. Using the above 
parameters in (|17|). we get z = 22% and 5 = 134° from PQCD calculation, which shows 
the dominance of the tree contribution in this decay and a large strong phase calculated 
from PQCD. 

Similarly, the decay amplitude for B s — > n~K + can be parameterized as 

A = V ub V: d T - V tb V td P = V ub V: d T[l + ze^+% (20) 

Therefore the averaged decay width for B Q S (B° S ) — > tt ± K t is 

r 2 A/f 3 

mm - tt^) = ^(|ai 2 /2+iai 2 /2) 

c 2 m 3 

= -T^r l \V: b V ud T\ 2 [l + 2z C osa C os5 + z 2 }. (21) 

It is a function of cos a cos 5. 

In Fig. EJ we plot the averaged branching ratio of the decay B° S (B° S ) — > 7r ± ii' =F with 
respect to the parameter a. Since the latest experiment constraint upon the CKM angle 
a; from Belle and BaBar is a around 100° [13j, we can arrive from Fig. |21 

6.2 x 10~ 6 < Br(B° s (B° s ) -> tt ± ^ t ) < 8.1 x 10" 6 , for 70° < a < 130°. (22) 




25 50 75 100 125 150 175 



a(degree) 



FIG. 2: The averaged branching ratio of B®(B S ) — > ir^K^ decay as a function of CKM angle a. 



i gives a similar branching ratios 
14 1 . In paper Q], Beneke eta/ 



Previous naive and generalized factorization approacl 
at 6 - 9 x lCr 6 with the form factor F Bs ^ K ~ 0.27 
also calculate this decay mode using QCD improved factorization approach (BBNS). It 
is based on naive factorization approach. The dominant contribution is still proportional 
to B s — > K form factor, which is introduced as an input parameter. In principal, the 
decay amplitude expand as series of a s and A/mg. But in practice, only the first order 
of a s corrections is calculated, including the so called non-factorizable contributions. The 
annihilation type contribution is power (A/mg) suppressed in BBNS approach. Therefore, 
the branching ratio predicted in QCD factorization and PQCD should not differ too much; 
but the CP violation in these two approaches will be different, since it depends on many 



non- leading order contributions (See below for discussion). In Ref. [15|, the branching ratio 
is about 10 x 10~ 6 , which is larger than our PQCD result and previous FA method jl^ . 
because their form factor F*--*(0) = 0.31 Q is larger than the previous factorization 
approach and our calculation below. 

The diagrams (a) and (b) in Fig. Q correspond to the B s — > K transition form factor 
F B "~* (q 2 = m 2 ~ 0), where q = Pi — P 2 is the momentum transfer. The sum of 
their amplitudes have been given by Eq. (jSJ), so we can use PQCD approach to compute 
this form factor. Our result is F Bs ^ K (0) = 0.27, if u b = 0.5; and F Bs ^ K (0) = 0.32, 
if Ub = 0.45. In our approach, this form factor is sensitive to the decay constant and 
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wave function of B s meson, where there is large uncertainty; but not sensitive to the K 
meson wave function. Eventually this form factor can be extracted from semi-leptonic 
experiments B s — > K~l + vi in the future. 

In our calculation, the only input parameters are wave functions, which stand for the 
non-perturbative contributions. Up to now, no exact solution is made for them. So 
the main uncertainty in PQCD approach comes from B s ,K,n wave functions. In this 
)aper, we choose the light cone wave functions which are obtained from QCD Sum Rules 
la,ll2|. For Ti meson, the distribution amplitude of light cone wave function should take 
asymptotic form if the energy scale // — > oo. But in our case, the scale is not more than 
5GeV, so we choose the corrected asymptotic form for twist 2 distribution amplitude <j>£, 
and other twist 3 distribution amplitudes derived using equation of motion by neglecting 



three particle wave functions [12(. These functions are listed in the Appendix, which are 



also used in decay mode B — ► Kir |7( and B — > nrr [8] etc. 

We also try to use the asymptotic form for 7r meson, for all the three distribution 
amplitudes (b&, 0^ and 4>^, since we have very poor knowledge about twist 3 distribution 
amplitudes |l8j]. The branching ratio of B s — ► tt + K~ is nearly unchanged (only 3%), 
because the branching ratio of B s — > n + K~ is mainly determined by the form factor 
F *~* (0) (see Fig. 1(a) and (b)) which is not dependent on ir wave function. However, the 
CP asymmetry changes from —28% to —13% by —54%, when a = 100°. This is because 
the direct CP asymmetry depend on the strong phase (see discussion below), which comes 
from non-factorizable and annihilation diagrams, where all three meson wave functions 
are involved. The CP asymmetry predicted here should be used with great care, since it 
depends on two much uncertainties. 

For heav y B and B s meson, its wave function is still under discussion using different 
approaches In this paper, we find the branching ratio of B®(B S ) — > n ± K T is sensitive 
to the wave function parameter Ub- For 0.45 < < 0.5, the resulted branching ratio will 
decrease from about 10 x 10 -6 to about 7 x 10~ 6 . When we set u>b = 0.45, our result is 
more closer to that of QCD factorization Q|. This sensitive dependence should be fixed 
by the B s — > K form factors from the semi-leptonic B s decays. Other uncertainties in 
our calculation include the next-to- leading order a s QCD corrections and higher twist 
contributions, which need more complicated calculations. 

From our calculation, we find that the dominant contribution comes from tree level 
diagrams (see Fig.l (a) and (b)) in this decay. If SU(3) symmetry is good, the branching 
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FIG. 3: The averaged branching ratio of B S (B S ) — ► ir K (K ) decay as a function of CKM angle 
a. 

ratio of B s — ► 7r + i^~ should be equal to that of B° — > 7r + 7r~. The experimental result of 
5° -> 7T+7T- is Br(B -> 7T+7T-) = (4.3+^ ± 0.5) x 10" 6 |2(|. The predicted branching 
ratio of B s — > 7rfT is about 1.7 times that of — > 7r + 7r~, where the difference comes 
mainly from SU(3) symmetry breaking: the decay constant fs s larger than fg and fx 
larger than f w . In the calculation, we also find that the electroweak-penguins contribution 
is negligibly small as 0.001% in branching ratio. 

For the experimental side, there is recent upper limit on the decay B® — > n + K~ [21], 

Br(B° s -> -K+K-) < 7.5 x 10~ 6 , (23) 

at 90% C.L. Our predicted result is consistent with this upper limit. 

For the decays of B S (B S ) — > 7i°K°(K ), the tree level contribution is suppressed due 
to the small Wilson coefficients C\ + Cz/S. Thus the penguin diagram contribution is 
comparable with the tree contribution. We study the averaged branching ratio of the 
decay B S (B S ) — > ir°K°(K ) as a function of a in Fig. |3J It is similar with Fig|21 We find 
that the branching ratio of B S (B S ) 7r°K° (K ®) is about 1.8 x 10~ 7 when a is near 100°, 
it is a little smaller than the result of Ref . 15( . 

In SM, the CKM phase angle is the origin of CP violation. Using Eqs. ()19|) and (|2Uj) . 
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the direct CP violation parameter can be derived as 

a cp = hrb i i j\2 = i , ; o- ( 24 ) 





2 _ 


.4 


— 2z sin a sin 5 




2 + 




2 1 + 22; cos a cos 5 + z 2 



It is approximately proportional to CKM angle sin a, strong phase sin 5 and the relative 
size z between penguin contribution and tree contribution. We show the direct CP vi- 
olation parameters as a function of CKM angle a in Fig. 0] From this figure one can 
see that the direct CP asymmetry parameter of B®(B®) — > ir^K^ and 7r ^ (i^ ) can be 
as large as —31% and —62% when a is near 75°. The larger direct CP asymmetry of 
B° S (B° S ) -> 7T°K°(K ) decay is mainly due to a larger z in -> tt°K°(K ) than in 

The direct CP asymmetry predicted in QCD factorization approach is quite differ- 
ent from our result, due to the different source of strong phases. In QCD factorization 
approach, the strong phase mainly comes from the perturbative charm quark loop dia- 
gram, which is a s suppressed ^j]. While the strong phase in PQCD comes mainly from 
non-factorizable and annihilation type diagrams. The sign of the direct CP asymmetry 
is different for these two approaches in B Q S (B° S ) — > ix^K^ decay, and the magnitude of 
CP asymmetry in QCD factorization (about 5%) is also smaller than PQCD. The future 
LHC-b experiments can make a test for the two methods. 

For the decays of B S (B S ) — > -k K°(K ), the final K°(K°) mesons can not be detected 
directly. What the experiments measured are their mixtures K s and Kl, thus a mixing 
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induced CP violation is involved. Following notations in the previous literature 
define the mixing induced CP violation parameter as 

-2/m(A cf ) 

a < + « = i + iw (25) 

where 

= Vj,V u {n K°\H eff \B°) 
CP V tb V t t(7r°K°\H eff \B°y 1 ] 

Using unitarity condition of the CKM matrix V tb V* d = -V ub V* d - V cb V c * d , and Eas. (ll9l2()l) . 

we can get 

\ CP = ±^L±, (27) 

e l T + x 

where x = ^r h X?i x +°g. Combining eq. (l2Tl) and (1231) . we can get 

I V ub v u d\ 1 +" 

sin 27 + 2Re(x) sin 7 



1 + \x\ 2 + 2Re(x) cos 7 



(28) 
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FIG. 4: Direct CP violation parameters of B®(B S ) — > ir^K^ (dashed line) and B S (B S ) — > 
7T K°(K ) (solid line) as a function of CKM angle a. 

If \x\ is a very small number, the mixing induced CP asymmetry is proportional to sin 27, 
which will be a good place for the CKM angle 7 measurement. However as we already 
mentioned, the tree contribution in this channel is suppressed, \x\ = 2.3 is a large number, 
so that the sin 7 behavior is dominant in the eq. (|28)1. The result of mixing induced CP 
violation is shown in Fig. which is indeed a roughly sin 7 behavior. The tail near 
7 ~ 180° also shows the contribution from sin 27 in eq. (J28j) . 

IV. SUMMARY 

In this work, we study the branching ratio and CP asymmetry of the decays B®(B®) — > 
■k^K^ and B S (B S ) — > tt K°(K ) in PQCD approach. From our calculation, we find 
that the branching ratio of B°(B° S ) tt ± K t is about (6 ~ 10) x 1(T 6 ; Br(B s (B s ) 
7T K°(K )) around 2 x 10~ 7 and there are large CP violation in the processes, which may 
be measured in the future LHC-b experiments and BTeV experiments at Fermilab. 
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FIG. 5: Mixing induced CP violation parameter of B S (B S ) — ► ir° K° (K ) as a function of CKM 
angle 7. 
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APPENDIX A: FORMULAS FOR THE CALCULATIONS USED IN THE 
TEXT 



In the appendix we present the explicit expressions of the formulas used in section II. 
First, we give the expressions of the meson distribution amplitudes 4>m- For B s meson 
wave function, we use the similar wave function as B meson 0, 0] : 



4>B s (x,b) = N Ba x 2 (l -x) 2 exp 



M 2 B x 2 



: (u b bf 



(Al) 



We set the central value of parameter uo b = 0.5 GeV in our numerical calculation, and 
Nb b = 63.7GeV is the normalization constant using fs a = 230MeV. 
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The 7r meson's distribution amplitudes are given by light cone QCD sum rules J 1 7[ : 



- x) {l + 0.44C 2 3/2 (t) + 0.25Q 3/2 (t)} 
J:: f l + 0.43C 2 1/2 (t) + 0.09Q 1/2 (t)} , 



2y/2N c 



1 - 2x) {1 + 0.55(10x 2 - lOx + 1)} 



2y/2W c 

where t = 1 — 2x. The Gegenbauer polynomials are defined by: 



(A2) 



1 



Cl /2 {t) = ^(35t 4 -30t 2 + 3), 



Cl ,2 {t) = ^(5t 2 - 1), Cl /2 {t) = ^(21t 4 - 14t 2 + 1). 



(A3) 



APT- 

We use the distribution amplitude <f) K ' ' of the K meson from Ref. 
6/x 



[16]: 



.r 



2y/2N~ c 

fx 



2^2N C 
fx 



x(l - x)[l + O.lbt + 0.405(5t 2 - 1)], 
[1 + 0.106(3t 2 - 1) - 0.148(3 - 30t 2 + 35t 4 )/8] 
t[l + 0.1581(5t 2 -3)], 



(A4) 



2V2N, 

whose coefficients correspond to itiqk = 1.6GeV. 

In our numerical analysis, we use the one loop expression for the strong running cou- 
pling constant, 



(A5) 



A)hi(/i 2 /A 2 ) 

where (3q = (33 — 2n/)/3 and n/ is the number of active quark flavor at the appropriate 
scale p. A is the QCD scale, which we take A = 250MeV at n/ = 4. 
Sb s , Sk- used in the decay amplitudes are defined as 

** dp 



S Bs (t) = s(x 1 P 1 + 1 b 1 ) + 2 [ ^(a s (p)), 

Ji/biP 

S n+ (t) = s(x 3 P+, b 3 ) + s((l - x 3 )P+, b 3 ) + 2 f 

Ji. 



(A6) 
(A7) 



S K -(t) = s(x 2 P^,b 2 )+s((l-x 2 )P 2 -,b 2 ) + 2 / ^(asifl)), (A8) 

where the so called Sudakov factor s(Q, b) resulting from the resummation of double 
logarithms is given as 0, 



^7(a s (/i)), 

l/b 3 ^ 

1 dp 



s(Q,b) 



hi 



A(a(p)) + B{a s (p)) 



(A9) 
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with 



A : 

B 





"67 


c F ^ + 


9 


7T 




e 2 7B -l 






3 7T V 


2 



7T" 

y 



10 



27 



n/ + ~/3 Q In — 



7T 



(A10) 
(All) 



Here 7^ = 0.57722- • • is the Euler constant, n/ is the active quark flavor number. For 



the detailed derivation of the Sudakov factors, see Ref. 

The functions hi(i = a,c.e.g) come from the Fourier transformation of propagators of 
virtual quark and gluon in the hard part calculations. They are given as 



h a (xi,X2,bi,b 2 ) = S t (x2)K (M B y/xix 2 b 1 ) 
x [9(b 2 - 6 1 )/o(M Bv ^6 1 )ir (M Bv ^6 2 ) + (&! 



(A12) 



h^\xx,x 2 ,X3,b 2 ,b 3 ) = 
0(b 2 - b 3 )I (M B y/ Xl (l - x 2 )b 3 )K (M BV / Xl (l - x 2 )b 2 ) 

K (M B F (j) b 3 ), for Fh > o' 
+ (63 <-> 63) )■ x I . m ? I , (A13) 



where Hq 1 (z) = Jo (z) + iY (z), and Fq^s are defined by 

-F( 2 !) = xx + x 2 + x 3 - xix 2 - x 2 x 3 - 1, F 2 2) = xi - x 3 - xix 2 + x 2 x 3 ; (A14) 

h<J\x 1 ,x 2 ,x 3 ,b h b 2 ) = 
0(b 2 - 6i)yHi 1) (M Bv ^6 2 )Jo(M Bv ^6 1 ) 

^^^H^^^for ^,<oJ' (A15) 
where -F e (j)'s are defined by 

F e(i) = xi + x 2 + x 3 - xix 2 - x 2 x 3 , F 2 (2) = xix 2 - x 2 x 3 ; (A16) 

h g (x 2 ,x 3 ,b 2 ,b 3 ) = S t {x 2 ) 1 ^-H^\M B ^x 2 x 3 b 3 ) 

x [0(63 - b 2 )UM B ^ 2 b 2 )^-H^{M By /r 2 b 3 ) + (6 2 <- 63)]- (A17) 
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We adopt the parametrization for S t {x) contributing to the factorizable diagrams 261 ] . 

01 + 2CTV0/O I \ 

St(x)= rr \ M \ V (l-*)] c , c=0.3. (A18) 
The hard scale t^s in Eq.([H])- (jl2)) are chosen as 

t\ = m^MB^T^^l/b^l/h), 

t\ = max (M Bv /x7, 1/6 1? l/6 2 ), 

*J = max(M Bl ^y , M BV /xi(l - x 2 ), l/6 2 , l/6 3 ), 

^ = max(M B ^F^, M B ^(1 - x 2 ), 1/&2, 1/63), 

t* = max(M B S J\F^\,M B ^/x^, l/b 1 ,l/b 2 ), 

t\ = max{M B ^\F^ 2) \, M B ^T^ 1/h, l/b 2 ), 
t\ = max(M B ^/x 2 ',l/b 2 ,l/b 3 ), 

t\ = max(M Bv ^, 1/62,1/63). (A19) 

They are given as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 
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